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Analysis  ^of  the  aerodynamics  and  jet  reaction  control 
of  a 5 "/54  gun-launched  missile  (GLM)  was  conducted.  Aero- 
dynamic analysis  included  estimation  of  list  and  drag  coef- 
ficients, lift-to-drag  ratio,  and  center  of  pressure  location 
by  calculation  for  several  versions  of  the  missile  at  Mach 
numbers  between  Mach  2.0  and  3. u.  Analysis  included  emphasis 
on  compatibility  with  existing  5"/54  Mk  45  gun  system.  Jet 
reaction  control  system  investigation  included  calculation 
of  the  side  force  required  to  trim  the  missile  for  a 30-g 
maneuver  at  Mach  3.0. 
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I . INTRODUCTION 


Since  the  attack  and  destruction  of  the  Israeli  destroyer 
"Ela.-'-h"  in  1968  by  four  Soviet  built  SSN-2  Styx  missiles  fired 
from  Egyptian  missile  boats,  the  U.S.  Navy  has  been  searching 
for  combat  capabilities  to  counter  the  anti -ship  cruise  mis- 
sile threat.  In  this  effort,  the  capabilities  of  the  corabau 
systems  developed  have  ranged  from  the  close-in  self  defense 
systems  with  kill  range  capabilities  of  only  several  miles 
to  the  recently  developed,  multi-faceted  multi-threat  Aegis 
system  with  its  SM-2  missiles  capable  of  intercepting  poten- 
tial missile  threats  60  miles  [11  from  the  firing  ship. 

In  this  twelve  year  effort,  one  area  the  Navy  has  failed 
to  exploit  fully  is  the  development,  production,  and  deploy- 
ment of  a combat  system  that  CU  has  the  capability  to  inter- 
cept and  destroy  the  cruise  missile  threat  with  a reasonable 

P (probability  of  kill)  in  the  current  and  near  term  envi- 
K 

ronment,  (2)  has  a wide  and  immediate  applicability  to  the 
currently  existing  fleet  assets,  (3)  has  a medium  self- 
defense  range  capability,  and  (4)  Still  has  affordability  vis 

a vis  the  existing  alternatives. 

A combat  capability  that  may  meet  the  four  criteria  men- 
tioned previously  is  the  5"/54  gun-launched  missile  (GLM) 
currently  being  investigated  by  offices  in  the  Naval  Sea 
Systems  Command,  this  student,  and  private  contractors. 
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The  5"/54  GLM  is  a relatively  low  cost  combat  system 
with  immediate  applicabilicy  to  e^ery  5"/54  Mk  45  capable 
ship  in  the  fleet  with  or  without  a surface  to  air  capa- 
bility.. The  GLM  will  increase  significantly  the  anti-ship 
missile  defense  ( ASMD)  capability  over  the  conventional 
anti-air  gunnery  abilities  currently  existing. 

The  purpose  of  this  thesis  is  to  investigate  two 
major  aspects  of  the  GLM.  The  first  major  section  addresses 
the  aerodynamics  of  several  configurations  proposed  as  gun- 
launched  missiles.  Refer  to  the  theses  of  Parks  [2],  Frazier 
[3],  and  Brown  [4]  for  additional  information  concerning  the 
GLM.  The  second  major  section  addresses  the  field  of  jet 
reaction  control  as  a potential  control  system  in  a GLM. 
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I.  BASIC  ASSUMPTIONS  AND  DESIGN  CONSTRAINTS 


Basic  aerodynamic  and  reaction  jet  control  assumptions 
and  decisions  for  the  5”/54  caliber  gun-launched  mis-^le  (GLM) 
were  made  at  the  beginning  of  this  thesis  investigation. 

These  assumptions  and  decisions  were  made  on  the  basis  of 
mission  applicability  and  subsystem  requirements  and  are 
described  below. 

Ground  Rules;  Ground  rules  were  defined  for  both  the  5"/54 
(GLM)  and  the  gun;  the  ground  rules  are  now  discussed. 

Missile ; The  missile  flight  profile  was  considered  to  be  as 
follows : 

. Launch,  with  immediate  de-spin  at  muzzle  exit 

. Powered  ballistic  flight  throughout  trajectory  to  inter- 
cept, utilizing  the  concept  of  thrust  equal  to  drag 

. Launch  velocity  at  the  speed  corresponding  to  no  change 
in  the  powder  case  propellant  characteristics,  i.e. 
approximately  2100  fps 

. Missile  maneuver  capability  of  30  g's  at  Mach  3 and 
sea  level 

. Flight  time  of  approximately  20  seconds  and  range  limits 
of  10  nm  on  the  trajectory 

. Reaction  jet  control  or  combination  tail/reaction  jet 
control 

. Maximum  angle  of  attack  limited  only  by  the  axially 
symmetric  inlet 

Gun : In  accordance  with  the  physical  characteristics  of  the 

5"/54  caliber,  Mk  45  lightweight  gun,  the  missile  was  to  be 
no  greater  than  61  inches  in  length  and  be  able  to  interface 
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With  existing  loading  system  constraints.  The  gun  constraints 
force  a bottleneck  configuration  in  the  forebody  or  the 

missile. 

Subsystem  Assumptions:  For  realistic  aerodynamic  «vd  jet 
control  configurations,  basic  assumptions  made  about 

the  types  and  characteristics  of  the  ncn-aero«ynamic  subsystems 
Some  of  the  subsystems  are  dealt  with  extensively  in  the 
companion  reports  [2,  3,  4]  to  this  thesis.  The  subsystems, 
however,  are  addressed  to  some  extent  below,  as  they  pertain 
tc  the  aerodynamic  and  control  packages  in  this  thesis. 
Propulsion:  Integral  rocket  ramjet  propulsion  system  was 

assumed  with  ignition  immediately  after  launch. 

Warhead:  A fragmenting  warhead  capable  of  defeating  an  air 

target  with  the  capabilities  predicted  for  the  1980-1990 
time -frame  was  assumed. 

Structure-.  An  overall  structural  ability  to  withstand  launch 
from  the  5"/54  gun  was  assumed. 

Sensor:  An  active  sensor  system  located  in  the  forebody  of 

the  missile  capable  of  opposing  an  air  target  with  the  capabil 
ities  predicted  for  the  1980-1990  time  frame  was  assuirt**- 
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HI.  AERODYNAMIC  THEORY 


i 


A.  INTRODUCTION  TO  AERODYNAMIC  DESIGN 

In  determining  the  optimum  aerodynamic  properties  of  a 
particular  flight  vehicle,  the  aerodynamicist  has  several 
options.  One  method  is  construction  of  scale  models  of  the 
various  proposed  designs  for  testing,  with  the  result  being 
a determination  of  the  optimum  configuration,  aerodynamically , 
to  meet  the  mission  requirements  and  still  interface  with  all 
other  system  requirements  and  assumptions.  Unfortunately, 
even  scale  models  have  become  increasingly  cost  prohibitive 
in  these  days  of  high  inflation;  and  at  the  very  least,  the 
model  construction  consumes  considerable  time.  More  frequent- 
ly, the  aerodynamicist  undertakes  a feasibility  study  analy- 
tically, using  existing  aerodynamic  and  gas  dynamic  theory, 
calculating  the  required  parameters  of  lift,  drag,  center  of 
pressure,  etc.,  based  on  the  particular  shapes  imagined  for 
the  prescribed  mission.  Then,  after  iteration  with  the  other 
requirements  of  the  flight  vehicle,  a configuration  design 
is  built  and  tested. 

In  this  thesis,  the  aerodynamic  configurations  of  several 
options  for  the  5"/54  gun-launched  missile  (GLM)  were  under- 
taken through  aerodynamic  theory. 

Calculations  for  each  particular  area  of  concern  for  the 
missile  were  conducted.  Finally,  the  calculated  values  for 
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each  proposed  GLM  configuration  were  presented. 

For  each.  GLM  configuration,  the  particular  vehicle  was 
divided  into  four  of  the  five  major  areas  which  affected  the 
drag  and  lift  of  the  vehicle.  The  five  major  areas  which 
affect  the  drag  and  lift  of  a given  flight  vehicle  are  1)  the 
nose,  2)  the  body,  3)  the  wings  or  fins,  4)  the  base,  and  5) 
aerodynamic  interference.  Each  component  contributes  to  the 
overall  drag  and  lift  characteristics  which  ultimately  deter- 
mine the  aerodynamic  capability  of  the  missile.  Aerodynamic 
interference,  which  is  discussed  by  Nielsen  [5],  has  not  been 
considered. 

In  this  initial  design  effort,  three  iterative  aerodynamic 
configurations  were  chosen  for  study  and  evaluation;  each 
configuration  net  the  general  assumptions  and  ground  rules 
discussed  ear  er  for  GLM  application.  The  configurations 
were  propelled  by  ramjet  engines  using  axially  symmetric  nose 
inlets.  ,v""  erodynamic  forces  were  divided  into  two  areas 
of  conce”  external  aerodynamics  and  internal  aerodynamics. 
This  the  "insiders  external  aerodynamics. 

B.  NOSE  L .T 

Recall  from  the  introduction  to  aerodynamic  theory  that 
the  nose  of  the  flight  vehicle  is  a contributor  to  the  aero- 
dynamic lift  of  a body. 

Since  the  nose  inlet  configurations  were  dictated  by  the 
ramjet  engine  propulsion  system,  the  lift  due  to  this  particu- 
lar geometry  had  to  be  considered.  Hoernei  [6]  states  that 
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lift  is  that  force  which  is  due  to  the  pressure  distribution 
over  the  body  in  a direction  normal  to  the  freestream  veloc- 
ity. 

For  both  GLM  Configurations  I and  II,  Figure  1 depicts 
the  chosen  inlet  geometry.  Although  the  overall  dimensions 
were  slightly  different  in  each  configuration,  the  relative 
configurations  are  the  same  for  aerodynamic  purposes. 


Figure  1.  Forebody  and  inlet  of  Configurations  I and  II. 

To  analyze  and  determine  the  lift  generated  by  the  inlet 
nose  of  the  GLM,  Newton's  second  law  of  motion  is  needed. 
Considering  motion  of  an  isolated  body,  the  law  states  that 
the  summation  of  the  forces  exerted  on  a body  at  a certain 
instant  is  equal  to  the  rate  of  ch  in  momentum  of  the 

body  at  that  instant.  Hoerner  [6  ver,  modifies  the 

situation  in  that  for  a stream  of  air,  there  is  not  an  iso- 
lated finite  mass;  rather,  there  is  a mass  flow.  Momentum 
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is  transferred,  as  one  would  expect,  onco  the  passing  stream 
of  air  by  adding  a vertical  component  of  velocity. 

In  order  to  determine  the  force  created  by  the  nose 
inlet,  the  inlet  was  modeled  analytically  as  a control 
volume  with  a constant  mass  flux  and  a change  in  direction 
of  flow  equal  to  the  angle  of  attack  a.  Figure  2 depict 


the  dynamics. 


Figure  2.  Analytical  model  (or  nose  lift. 

From  the  geometry  of  Figure  2 and  the  momentum  equation, 
the  lift  L is 

L = wVsina 

where  w is  the  mass  flow  rate,  kg/sec:  and  V is  the  velocity 
of  the  air,  m/sec.  For  very  small  angles,  the  sine  of  the 
angle  may  be  approximated  by  the  angle.  For  small  angles 
of  attack  o,  si no  is  approximately  a.  Therefore 
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By  the  continuity  equation1,  the  mass  flow  rate  w is  equal  to 

3 


the  product  of  the  density  p of  fluid,  kg/m  , the  area  S 
through  «hioh  the  fluid  passes,  m2 , and  the  velocity  V of  the 
fluid.  By  inserting  the  values  of  the  product  equal  to  w 

into  equation  (1)  , 

L = pV2Sa 

Hoerner  [6]  states  that  lift  can  be  represented  as  the  product 
of  the  dynamic  freestream  pressure,  q,  N/m2 , the  area  S 
previously  defined  and  a non-dimensional  term  called  the 
coefficient  of  lift,  CL-  The  dynamic  freestream  pressure  is 

shown  in  equation  (4) . 


L = CL  qS 


(3) 

q = >lPV2  (4) 

By  substituting  equation  (4)  into  equation  (3),  the  general 


term  for  the  coefficient  of  lift  can  be  determined.  This 
is  shown  in  equation  (5) . By  substituting  equation  (2)  into 
equation  (5) , the  coefficient  of  lift  can  be  determined  for 
the  nose  section  of  the  GLM. 


CL  = 


pV2S/2 


= ,_PV  Sa — = 2a 
L PV2S/2 


(5) 

(6) 


XFor  more  information  on  the  continuity  equation, 
Liepmann  and  Roshko  [7]  is  an  excellent  reference. 
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For  the  configurations  studied  in  this  thesis,  CL  for  a 
given  angle  of  attack  a are  tabulated  in  Section  V. 

For  Configuration  III,  Figure  3 depicts  the  chosen  inlet 
geometry.  The  aerodynamic  cowling  was  a result  of  the  design 
process.  Not  all  of  the  incoming  volume  of  air  could  be 
utilized  by  the  ramjet  engine;  consequently,  the  capture 
streamtube  area,  AQ,  was  reduced.  For  the  iteration  from 
Configuration  I to  III,  the  dimensions  were  sufficiently 
different  to  warrant  investigation  of  alternative  methods 
for  calculating  the  lift  coefficient  resulting  from  che  nose 
configuration.  The  major  geometrical  difference  between 
Configuration  I and  Configuration  III  is  the  cowling  shown 
in  Figure  3.  The  lift  contribution  due  to  the  air  which 
enters  the  ramjet  was  calculated  using  equation  (6) . For 
the  calculation  of  the  cowling  lift,  the  computer  program  of 
Tillotson  and  Schonberger  [8]  was  used.  For  an  angle  of 
attack  of  5 degrees,  the  program  gave  CL  due  to  the  cowling 
significantly  less  than  the  other  contributions  to  CL. 


Figure  3.  Nose  and  forebody  of  Configuration  III  with 
redesigned  cowling. 
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Consequently,  the  cowling  lift  was  ignored  in  subsequent  cal- 
culations. AS  a word  of  caution,  the  importance  of  cowling, 
lift  rhould  be  investigated  in  more  detail.  The  lift  coeff 
cient  for  given  angles  of  attack  also  is  tabulated  in  rhe 

design  section. 

C . NOSE  DRAG 

For  any  aerodynamic  configuration,  the  nose  or  forebody 
can  cause  significant  drag.  Since,  in  this  design  study, 
the  GLM  configurations  were  powered  by  ramjet  engines  with 
an  axially  symmetric  inlet,  the  inlet  geometry  was  chosen 

such  that  the  nose  drag  was  minimal,  i.e.,  the  mass  flow 

_ . , j • n rr\iT“©s  1 and  3 was  as  closs  to 
ratio  A0/Ac  as  depicted  in  Figures 

unity  as  possible.  Thus,  for  both  Configuration  X and  II, 
the  drag  coefficient  for  the  nose  was  assumed  to  be  equal  to 

zero  at  an  angle  of  attack  a equal  to  zero. 

After  the  initial  design  efforts,  however,  modification 
of  the  nose  design  was  necessary  as  the  mass  flow  ratro  A0/Ac 
depicted  in  Figure  1 was  not  suitable  for  the  ramjet  engine. 
Figure  3 depicted  the  revised  design.  Calculation  of  the  nose 
drag  due  to  the  cowling  was  made  based  on  a computer  program 
developed  by  Tillotson  and  Schonberger  [3]. 


D . BODY  DRAG 

The  body  drag  of  a flight  vehicle  is 
major  factors:  lrag  due  to  the  friction 
drag  as  a result  of  the  configuration  or 


composed  of  two 
of  the  body  and 
shape  of  the  body. 
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in  this  thesis , Configuration  I total  drag  was  affected  only 
by  the  friction  drag.  Configurations  II  and  III  were  affected 
by  both  factors  of  body  drag.  Both  form  drag  (pressure)  and 
skin  friction  drag  were  considered  in  the  aerodynamic  analysis. 

Hoerner  [9]  states  that  the  pressure  drag  is  that  force 
which  represents  the  streamwise  component  of  the  pressure 
integrated  over  the  entire  configuration.  In  contrast,  the 
skin  friction  drag  is  that  force  obtained  by  integrating  the 
streamwise  component  of  the  shear  stress  over  the  flight 
vehicle.  in  order  to  calculate  the  skin  friction  drag  of  a 
particular  body,  the  following  equation  is  defined: 

Drag  = /lnyds 

where  T is  the  local  skin  friction,  N/m2,  ny  is  the  local 
y-compcnent  of  the  normal  vector  to  the  surface  and  ds  is  an 
element  of  area,  m2.  Fluid  flow  is  in  the  x-direction.  In 
this  thesis,  an  approximation  of  the  equation  for  skin  fnc- 
tion  was  made. 

Drag  = xSr 

In  equation  (7) , Sr  represents  the  reference  area.  The 
approximation  was  a result  of  averaging  the  shear  stress 
over  the  entire  surface  of  the  body.  The  average  shear 
stress  t is  defined 

T _ „ r (8) 

1 = q Cfc 

where  q i/  the  dynamic  freestream  pressure  previously  defined 
on  page  19  and  Cfc,  the  coefficient  of  friction  for  compress- 
ible fluid  flow.  For  more  information  on  the  momentum  ard 
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energy  equations  that  define  the  dynamic  pressure,  the 
reader  is  referred  to  references  [7,  10,11]. 

Recall  that  for  a perfect  gas 

p = pRT 

and  the  speed  of  sound  is  defined  as 


where 


a = yRT  (10) 

y = ratio  of  specific  heats 

m2 

R = specific  gas  constant,  ■ 


Since  the  definition  of  Mach  number  is 

M = V/a  (11) 

by  substituting  equations  (9) , (10)  and  (11)  into  equation 

(4),  the  dynamic  pressure  becomes 

q = ypM2/2  d2) 

In  order  to  calculate  the  average  shear  stress  over  the  body 
and  ultimately  the  skin  friction  drag,  the  aerodynamicist 
must  determine  Cfc-  Usually  this  coefficient  is  determined 
by  finding  the  coefficient  of  friction  in  an  incompressible 
flow  and  then  adjusting  for  compressible  flow.  In  discussing 
fluids  of  both  an  incompressible  and  compressible  nature, 
viscosity  should  be  addressed.  Shear  stress  can  be  defined  as 
proportional  to  the  rate  of  angular  deformation  of  a given 
fluid.  The  proportionality  constant  in  the  equation  is  defined 
as  the  viscosity  of  the  fluid.  Dimensions  for  viscosity  are 

Newton-seconds  per  square  meter. 

The  ratio  of  inertia  forces  to  viscous  forces  on  a body_ 

in  a fluid  flow  is  defined  as  the  Reynolds  number 
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Re 


(13) 


P V L 

y 

L = the  appropriate  characteristic 
length  of  the  vehicle  under 
analysis,  m 

y = viscosity  of  the  fluid  through 
which  the  vehicle  is  passing, 

N - S/m2 


Since  the  fluid  through  which  the  GLM  passes  is  air,  the 
viscosity  of  air  for  the  particular  altitude  or  altitudes 
at  which  the  GLM  operates  is  needed.  Figure  4 depicts  the 
Reynolds  number  per  characteristic  length  L as  a function  of 
velocity  and  altitude  for  U.S.  Standard  Atmosphere.  The 


characteristic  length  L was  equal  to  one  meter 


Figure 


4. 


Reynolds 

function 

Standard 


number  per  characteristic  length  as  a 
of  velocity  and  altitude  for  U.S. 
Atmosphere.  From  Bertin  and  Smith  (10 J 
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Hoerner  [9]  states  that  outside  of  the  boundary  layer, 
the  effects  of  the  shear  forces  are  negligible.  Because  of 
this,  t.ie  region  outside  of  the  boundary  layer  is  treated  as 
inviscid  flow.  For  airfoils  passing  through  air,  a complex 
flow  pattern  develops  as  the  velocity  of  the  flight  vehicle 
increases.  The  boundary  layer  is  divided  into  two  portions: 

(1)  a laminar  portion  and  (2)  turbulent  portion.  In  general, 
the  laminar  portion  of  the  boundary  layer  is  that  region  where 
the  "layers"  of  gas  do  not  mix  with  each  other  and  are  parallel 
to  each  other  with  steady  velocities. 

The  turbulent  portion,  in  contrast,  as  Hoerner  [9]  states, 
is  that  portion  in  "a  constant  'state  of  commotion  and  agita- 
tion', consisting  of  velocity  fluctuations  superimposed  to 
the  main  flow.  Turbulent  flow  occurs  within  boundary  layers 
(at  higher  Reynolds  numbers)  and  within  the  wake  behind  solid 
bodies."  For  vehicles  at  high  Reynolds  numbers  and  high  velo- 
cities, the  flow  is  characterized  by  a laminar  portion  of  the 
boundary  layer  that  then  transitions  through  a finite  region 
to  a turbulent  portion  of  boundary  layer.  Figure  5 illustrates 
this  complex  transition. 

T <YJ). 

•itlra 


X. 

Eigure  5.  Transition  from  laminar  to  turbulent  fluid  flow  over 
airfoil.  From  Talay  [12] 
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(14) 


Thus,  in  the  laminar  region 


Di  ~ Cfc  qSl 


and  in  the  turbulent  region 


Dt  = CftqSt 


(15) 


The  coefficient  of  friction  for  laminar  incompressible  flow 

has  been  determined  empirically  [6,  9]  to  be 

C = 1.328/  /Re"  {16) 

in  a similar  manner,  Schoenherr,  in  Hoerner  [9],  established  an 
equation  which  is  used  for  calculating  the  urfculent  coeffi- 
cient of  friction  for  incompressible  flow  and  i ip  to 

the  highest  Reynolds  numbers  normally  encounter*  d in  .gineer- 
ing,  i.e.,  around  109.  The  Schoenherr  equation  is 

log10(ReCf)  = 0.242/C^  <17) 

For  a given  Re,  an  iterative  process  is  conducted  until  both 
sides  of  the  equation  are  balanced,  which  determines  the  Cf. 

Thus,  the  drag  due  to  skin  friction  on  the  body  of  the 
flight  vehicle  in  incompressible  flow  is 


D — qC^  St  + qC^ 


(18) 


If  x*  is  defined  as  the  distance  from  the  leading  edge  or 
nose  of  the  body  to  the  point  at  which  laminar  flow  changes 
to  turbulent  flow,  and  if  the  surface  area  of  the  entire  body 

is  S (L)  , then 

Drag  = q CfiS(x*)  + qCft  [A(L)-A(x*)]  (19) 

since  A(L)  - A(x*)  would  be  the  surface  area  of  the  turbulent 

region . 


The  distance  x*  along  the  body  at  which  the  transition 
occurs  is  given  by 

* U Re*  (20) 

X " p V 

where  Re*  is  the  Reynolds  number  at  transition. 

For  Reynolds  numbers  greater  than  1 or  2 million,  based 
on  the  total  length  of  a given  body,  Hoerner  [9]  plots  a Re* 
of  2 x 10^. 

Once  the  coefficient  of  friction  for  an  incompressible 
flow  is  found,  it  is  adjusted  " r compressibility  of  air, 
since  the  GLM  is  traveling  ically  through  the  air. 

An  equation  from  Hoerner  [9]  * jsed  for  adjusting  from 

incompressible  to  compressible  flow. 

c,  /Cf.  = (1  ♦ O.JSM2)-0-58  (21) 

f C tl 

Calculation  of  a coefficient  of  friction  for  compressible 
flow  at  a given  Mach  number  is  accomplished  by  using  the 
discussed  earlier. 

Finally,  returning  to  equation  (8) , 

T " * Cfc 

the  average  shear  stress  over  the  surface  of  the  body  is 
calculated;  and  using  equation  (7) , the  drag  of  the  body 

due  to  friction  is  determined. 

For  configuration  I,  in  which  the  body  drag  was  due 
solely  to  skin  friction,  the  calculations  are  discussed  in 
detail.  SI  units  were  used;  and  the  density,  pressure,  and 
temperature  values  for  the  GLM  at  sea  level  [6]  were  used  in 
the  calculations. 


27 


Since  several  configurations  were  investigated  in  this 
thesis,  and  each  was  evaluated  at  various  velocities  and 
altitudes , a computer  program  using  the  HP  9830  was  developed. 
This  program  is  discussed  in  greater  detail  in  Section  I/. 

E.  BODY  LIFT 

The  missile  designer  chooses  a slender,  streamline  shape 
to  minimize  the  drag  on  the  body,  as  has  been  shown  already 
in  the  section  on  body  drag.  One  additional,  positive  side 
effect  is  that  the  streamline  shape  also  produces  a lift  on 
the  body.  Body  lift  was  calculated  for  the  GLM  configurations 
investigated.  According  to  Hoerner  [6],  for  streamline  bodies 
such  as  GLM  configurations  analyzed  in  this  thesis,  the  lift 
consists  of  two  separate  components,  one  due  to  circulation 

and  the  other  due  to  cross  flow. 

In  a slender  body  of  revolution,  such  as  the  cylindrical 
mid  section  of  the  GLM  body,  the  lift  coefficient  for  low 

angles  of  attack  a is 

CT  = sin  (2a)  = 2sinacosa 

1j 

For  small  angles,  CL  re  uces  to 

Cj.  = 2a 

As  the  angle  of  attack  increases  above  a few  degrees,  the 
component  of  flow  normal  tj  the  body  becomes  more  signifi- 
cant. The  second  component  of  body  lift  is  proportional  to 
the  cross-flow  coefficient  known  as 

Cc  = Cn  = C90 
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for  several  of  its  names.  The  proportionality  will  be  shown 
explicitly  in  equation  (22).  Hoerner  [6]  states  that  at  very 
high  Reynolds  numbers,  such  as  Re  * 10^,  experimental  results 
indicate  CgQ  between  1.1  and  1.2  for  the  type  cylinder  depic- 
ted in  Figure  6.  The  data  for  the  results  are  presumed  to  be  for 
subsonic  flow. 


Figure  6.  Cross-flow  coefficient  versus  angle 
of  attack  for  cylinders  in  subsonic 
flow.  From  Hoerner  [6]. 


The  equation  used  to  calculate  the  normal  force  coefficient 
is 


CN  2a  + C90  Sr 


a 


(22) 


The  ratio  S^/S^.  is  the  ratio  of  the  presented  area  to  the 
reference  area.  For  more  information  on  the  theory  and  the 
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experimental  results,  the  reader  is  referred  to  Hoerner  [6]. 

As  shown  in  Appendix  C,  lift  L is  nearly  equal  to  the 
normal  force  N when  the  angle  of  attack  is  small  and  the  lift— 
to-drag  ratio  is  large.  Since  all  the  configurations  were 
considered  to  be  operating  in  the  low  angle  of  attack  regime 
with  large  L/D,  the  body  normal  force  was  considered  approxi- 
mately equal  to  the  body  lift  coefficient. 

Lift  coefficients  as  a result  of  the  body  normal  forces 
for  all  configurations  are  tabulated  in  the  design  section. 

F.  TAILFIN  DRAG 

The  drag  of  the  airfoil  is  a result  of  two  force  compo- 
nents. These  two  forces  are  called  wave  (or  form)  Irag  and 
-fiction  drag. 

The  theory  defining  friction  drag  has  been  discussed 
already.  Detailed  calculations  relevant  to  the  particular 
configurations  are  tabulated  in  Section  V.  Wave  drag  is  due 
to  a pressure  exerting  a drag  on  the  GLM  fins  and  therefore 

is  discussed  in  its  entirety. 

For  the  GLM,  a slender  double-wedge  airfoil  in  supersonic 
flow,  as  shown  in  Figure  7,  was  considered  since  the  flight 
range  was  Mach  2 . 0 to  Mach  3.0. 
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Figure  7.  Double-wedge  airfoil  at  angle  of  attack 


For  the  double-wedge  airfoil  which  is  depicted  in  Figure  7 
the  supersonic  flow  causes  the  formation  of  a shock  wave  that 
is  attached  to  the  leading  edge  of  the  airfoil.  In  addition, 
the  wave  is  planar,  and  the  downstream  flow  is  isentropic. 

A technique  called  the  shock  expansion  theory  was  used  [1G] 
to  calculate  the  section  drag  coefficient  and  similarly  the 
section  lift  coefficient  for  the  double-wedge  airfoil. 

Liepmann  and  Roshko  [7]  state  shock-expansion  theory  uses 
the  principle  that  the  oblique  shock  wave  and  simple  lsentro- 
pic  wave  furnish  the  building  blocks  for  analyzing  problems 
in  two-dimensional  supersonic  flow.  Key  to  the  theory  is 
attachment  of  the  shocks.  The  two-dimensional  theory  ignores 
wing  tip  unloading,  and  this  thesis  will  not  address  wing  tip 
unloading  for  the  GLM . More  information  on  shock-expansion 
theory  is  provided  in  Refs.  [7,  10,  11]  • 
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The  geometry  of  the  chosen  fin  for  the  GLM  was  repre- 
sented in  Figure  7.  The  calculations,  shown  for  example 
purposes  to  illustrate  the  theory,  used  the  geometry  of  this 
fin.  Sample  calculations  are  for  an  angle  of  attack  a equaj. 
to  5 degrees. 

From  Figure  7,  the  surface  in  region  1 is  parallel  to 
the  freestream  condition.  The  flow  does  not  turn;  and 
consequently,  the  flow  properties  in  region  1 are  the  same 
as  freestream.  The  pressure  coefficient  on  the  airfoil  sur- 
face bounding  region  1 is  zero.  Thus,  = 3.0,  the 

Prandtl-Meyer  angle  = 49.75°,  8^  = 0°  and  Cpl  = 0. 

Since  the  surface  of  the  airfoil  in  region  2 turns  away 
from  the  flow  in  region  1,  the  flow  accelerates  isen tropically 
in  going  from  region  1 to  region  2.  From  Figure  7 dv  = -d0; 
and  since  02  = -10°,  \>2  = ^ - A0,  and  v2  = 59.75°.  Using 
compressible  air  flow  tables  [14],  M2  = 3.58. 

To  calculate  the  coefficient  of  pressure  in  region  2,  the 
Prandtl-Meyer  relation  discussed  in  Liepmann  and  Roshko  [7] 
is  used.  Equation  (23)  gives  the  coefficient  of  pressure 

desired. 

C „ 111!?. = -L_  ( - 1)  (23) 

P2  tt/2)P.M„2  Y P« 

Since  the  flow  over  the  upper  surface  of  the  airfoil  is 

isen tropic,  the  stagnation  pressure  is  constant,  and  ptoo  - Ptl 

= pt2-  Therefore  equation  (23)  becomes 
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(24) 


2 


( 


P2  Pf 


- 1) 


(24) 


C 


p2 


Pt2  p°° 


Using  the  values  for  the  airfoil  considered,  Cp2  =-  0.0904. 

Turning  to  tbe  lower  surface  of  the  airfoil,  the  pressure 
coefficients  for  region  1 ani  4 may  be  calculated. 

When  supersonic  flow  encounters  a change  in  direction 
resulting  in  compression,  the  flow  decelerates,  and  a shock 
wave  occurs  [7,  10,11  ].  The  shock  wave  is  a discontinuity  in 
the  continuation  of  the  flow.  These  discontinuities  may  be 
large  or  small,  hence  the  notion  of  "strong"  or  "weak"  shock 
waves.  Since  the  shock  wave  is  usually  thought  of  as  a large 
pressure  wave, [11]  the  term  "pressure  jump"  is  used  as  a 

measure  of  shock  strength. 

Various  parameters  may  be  used  to  measure  the  shock 
strength.  Since  the  flow  through  the  shock  wave  is  adiabatic, 
[7,  10,  11],  the  entropy  must  increase  as  the  flow  passes 
through  the  shock  wave.  Two  of  the  parameters  by  which  the 
shock  wave  may  be  measured  are  the  shock  wave  angle  6 and  the 
deflection  angle  5.  The  following  relation  exists  between 
these  two  angles: 

(Y  + 1)  M.2 

cot  <5  = tan  3 [ 5 5 - 1]  (25) 

2 (M^  sin  p-1) 

Figure  8 depicts  the  appropriate  angles  and  symbols. 
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Figure  8.  Shock,  wave  angle  and  deflection  angle  in 
supersonic  flow. 


Since  the  deflection  angle  5 is  known  for  the  airfoil  under 
consideration,  the  shock  wave  angle  3 can  be  calculated. 

Using  equation  (25),  3 = 27.5°.  The  pressure  jump  can  be 
found  in  region  3 by  M^sinS  = Mn  and  the  tables.  For  the 
example  airfoil,  = 1.385,  and  P^P^  = 2.0  71. 

Using  equation  (24),  Cp3  = 0.1700. 

The  calculation  of  M3  can  be  accomplished  by  using 

equation  (26)  with  5 = 10  . 

= (y  -1)  M0O2sin2 3 + 2 (26) 

"3  [2y.M  2sin2  3-(y-l)]  sin2  (3-5) 

00 

Inserting  the  necessary  values  yields  M3  = 2.48;  and  from  the 
tables  v3  = 38.75°. 

Since  from  Figure  7,  A0  = 10°  between  region  3 and  region 
4,  v4  = 48.75°. 

Returning  to  the  tables  with  the  Prandtl-Meyer  angle, 

M = 2.95;  and  by  using  equation  (24)  once  more,  Cp4  = 0.0011 
4 
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Having  determined  the  pressures  on  the  individual  regions 
of  the  thin  double-wedge  airfoil,  the  section  drag  coefficient 
can  be  determined.  This  is  the  drag  coefficient  for  an  a 
foil  with  unit  span.  Using  the  fact  that  the  net  force  in 
any  direction  due  to  constant  pressure  acting  on  a closed 
surface  is  zero,  [10]  the  pressure  coefficients  may  be  summed 
around  the  surface  of  the  airfoil  as  in  equation  (27)  to  find 
the  section  drag  coefficient 


J — j.Y  C sin  8 

'OSO  ' P 


•a-  - JcoiTzL'- p “ <27> 

In  this  case  ca,  = 0.0075.  Once  the  section  drag  coefficient 
is  determined,  the  drag  coefficient  can  be  calculated  by 
multiplying  the  value  by  the  span  b,  measured  in  meters,  of  the 

airfoil. 


(28) 


For  the  GLM  fin  under  consideration,  Ca  - 0.0013. 

The  coefficient  of  drag  for  each  tailfin  is  based  on  the 
fin  area  as  the  reference  area.  For  this  value  to  be  applica- 
ble to  the  GLM  under  consideration,  it  must  be  converted  to  a 
drag  coefficient  based  upon  the  reference  area  for  the  GLM. 

As  discussed  earlier,  the  reference  area  fcr  the  GLM  is  the 
base  area.  Thus,  the  drag  coefficient  as  a resuit  of  the  sum 
of  all  four  tailfins  was  multiplied  by  the  tailfin  reference 
area  divided  by  the  base  area.  Results  for  the  tailfin  configu- 

tabulated  in  Section  V. 


rations  selected  are 


G.  TAILFIN  LIFT 


For  thin  airfoils,  such,  as  the  double-wedge  airfoil  to 
be  used  in  the  GLM  configuration,  the  section  lift  coefficient 
is  found  as  a result  of  the  same  technique  as  used  in  calculat- 
ing the  section  drag  coefficient.  Once  the  missile  designer 
determines  the  pressures  of  the  individual  regions  of  the  air- 
foil, the  lift  coefficient  for  the  section  of  the  airfoil 
shown  in  Figure  7 can  be  determined.  Equation  (29)  shows  how 
the  section  lift  coefficient  is  calculated. 

CT  . * T Y C cos 9 (29) 

L'  2 cos  6 p 

Once  the  section  lift  coefficient  is  determined,  the  lift 
coefficient  can  be  calculated  by  multiplying  the  section  lift 
coefficient  by  the  span  of  the  airfoil  under  consideration. 

As  was  discussed  in  the  section  on  tailfin  drag,  the 
coefficient  of  lift  for  each  tailfin  is  based  on  the  reference 
area  for  the  fin.  This  coefficient  also  had  to  be  converted 
to  a lift  coefficient  based  on  the  reference  area  of  the  GLM. 
Thus,  the  lift  coefficient  based  upon  the  lifting  surfaces  of 
two  tailfins  was  multiplied  by  the  factor  of  the  tailfin 
reference  area  divided  by  the  base  area.  Results  for  the 
tailfin  configurations  are  tabulated  in  Section  V. 


IV.  AERODYNAMIC  DESIGN  CALCULATIONS 


A.  INTRODUCTION 

This  section  of  the  thesis  Consider^  the  detailed  calcu- 
lations for  the  aerodynamic  values  for  one  of  the  several 
configurations.  Sample  calculations  are  illustrated  for  the 
representative  sections  contributing  to  lift  or  drag  on 
Configuration  I.  Additionally,  the  several  computer  programs 
utilized  in  this  thesis  to  aid  in  more  rapid  calculations  of 
subsections  of  each  configuration  are  discussed  and  presented. 

In  the  course  of  the  aerodynamic  analysis,  calculations 
for  several  sections  of  each  configuration  we,  4>  identical 
except  for  the  basic  parameter  involved.  The  calculations 
presented  in  this  section  are  designed  to  illustrate  and 
amplify  the  aerodynamic  theory  formally  presented  in  Section 
III.  Sample  calculations  presented  in  this  section  are 
representative  in  nature,  and  not  every  value  tabulated  in 
Section  V is  treated.  Rather  the  purpose  of  this  section  i* 
to  illustrate  explicitly  how  the  design  of  the  several  confi- 
gurations for  this  thesis  was  approached.  Methodology  is 
emphasized,  and  the  fact  that  in  all  analytical  endeavors 
some  plain  "number  crunching"  must  be  done  is  illustrated. 

B.  SECTION  CALCULATIONS 

Calculations  for  each  component  of  Configuration  I are 
presented  in  this  section.  Table  I provides  the  geometric 
data  for  Configuration  I.  Tab lei  I provides  the  standard 

atmospheric  values  used  in  many  of  the  calculations. 
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TABLE 

I 

CONFIGURATION 

I GEOMETRIC 

DATA  FOR  CALCULATIONS 

BODY: 

TAILFIN : 

Length,  overall,  ro 

1.3716 

Area  (2  panels),  m2  0.0121 

Diameter,  m 

0.1270 

Chord,  m 0.0635 

Volume,  m3 

0.016  4 

Span,  m 0.09525 

2 

Wetted  Area,  m 

0.5168 

Thickness  ratio  0.0875 

2 

Base  area,  m 

0.0127 

Aspect  ratio  1.499 

TABLE 

II 

STANDARD  Ar;. 

:mospheric  values  for  air  at  sea  level 

5 = 1.225  kg/m3 

T = 288. 2°K 

P = 101,300  N/m2 


\ 


1.4 


1.  Calculation  of  Body  Drag 

From  the  data  in  Table  I,  the  surface  area  of  the 

body  was  calculated, 

A = dL  = 0.5168  m*  (30) 

and  the  dynamic  freestream  pressure  was  determined 

q = !,  pv2  = 637,245  N/m2  (31! 

The  Reynolds  number  which  was  defined  in  equation  (Hi  , was 

found  to  be 

(32) 


( i .225)  (1020)  (1.29541 
7-5 


= 89.9  x 10 

Re  = — 


1.8  x 10 

Since  the  Reynolds  nu^er  is  greater  than  the  transrtion 
Reynolds  number,  the  body  has  turbulent  flow;  refer  to 
Section  III  for  a discussion  of  boundary  layer  transition 
Recall  from  the  theory  that  for  incompressible  flow  in  a 
turbulent  boundary  layer 


log  10  ( Re  Cf)  = 0 . 242/Cf 


33) 


where  C£  is  the  coefficient  of  friction.  Equation  (33'  is 
solved  iteratively  using  the  Reynolds  number  from  equation  (32) 

c = 0.00211 

For  the  adjustment  from  incompressible  to  compressible 

flow,  equation  (34)  is  used  where  M=  3.0. 

2 ” (34) 

C£c  - C£i  (1  t 0.15  M ) 

After  inserting  the  incompressible  value  in  equation  (34i, 
c was  determined  to  have  a value  of  0.001285.  Recall  from 
the  discussion  of  body  drag  that  the  shear  stress  X was 
defined  as  x = q C£c-  Inserting  numerical  values  gives 

T = 818.26. 
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423.2  N 


Finally,  using  equation  (7)  , drag  - 

and  CD  = D/qsref  = 0.0524. 

Modifying  a program  developed  for  the  NPS  course  AE  4705, 

[13]  the  above  steps  were  programed  for  the  HP  9830  computer. 
The  program,  is  listed  in  Table  III  with  a sample  output  in 
Table  XV.  With  inputs  of  velocity,  diameter,  length  and 
altitude,  the  program  calculates  the  body  skin  friction  drag 
of  a cylindrical  body. 
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TABLE  III 

HP  9830  PROGRAM  TO  CALCULATE  BODY  DRAG;  PROGRAM  STEPS 


I 


It*  F'F'INT  TH I :■  F'R  LRRM  lALUJLFiTE':'  FFiijE'. i_  d'UDY  DRUG  . 

20  PRINT  “THIS  PROG  Phtl  IR  STUPED  IN  FILE  0. 

30  PRINT  "" 

40  DEG 
50  D 1 = 1 O0 

60  D I P “PRINT  I NhUT'.'Y£'E  = 1 NO-0'  * 

70  INPUT  1 
80  PRINT 

90  PRINT  •FOR  iUPEPSOuIC  '•'ELOCITY  INSERT  0. 

100  PRINT 

110  DISP  “SEE  PRINTED  STATEMENT  CONCERNING  ' Ei_0'..i  fV  ? 

120  INPUT  K 9 
130  Dl”1 00 
140  D=0 . 12593 
150  P=D  2 
160  H=PI-RT2 
170  Hi =340 
ISO  7=1020 
190  M = V Hi 
200  H=7620 
210  Y = 0 

220  PRINT  "*--■* -r  * -r + + + + + 

230  PRINT  "ALTITUDE  3" 5 V" METER" 

240  P0= 1.225 
250  Z=R0*EXR‘  -V  - HI- 
260  N7  = 1 . 3E-05 
2 7u  P 1 = 2 E + 0 6 

2 3 0 X = < P:  1*M7 * ESP  V H R 0 * V 

290  L= 1.5494 

292  IF  K9=0  THEN  310 

3 1 0 3 2 = 2 - F'  I *L  + F! 

320  IF  L>X  THEN  350 

330  3 =33 

340  GOTO  420 

350  5 = 2*PI*R-«-X 

420  C 1 = 1 .323  '-SOP'  PI > 

425  C 1 =C  1 < ( 1 +0 . 45*  MT2  ■ r 0 . 25  > 

430  R2=R  1 * (L •••“:> 

440  PRINT  ” 

450  0 1 = C 1 
460  02=LGT (Ql*R2  1 
470  03=0212 
4 g 0 Q 4=0.05 8 5 6 4 0 3 
490  05 = FIBS  < 04-01  > 

500  IF  05<1 E-06  THEN  535 
510  IF  05(1 E-05  THEN  570 
520  IF  05(1 E-04  THEN  550 
530  0 1 =0 1 *0 1 Dl 
54U  GOTO  460 


I 


1 7 

Aj 


A 


550  r i -r t L r l ' j-^Dl  > 

St. Li  GOTO  460 

570  0 1 =0 1 +C  1 40-D  1 ' 

580  GOTO  460 

585  01=01  ■ . 1 +0.  15~N  r-2  • ru.  53  • 

580  i*|8=Rij- 1 t 2 * 0 . 5 -£,  F -V  H 1 
600  D7  = Q9*  •'  C 1 - :•  ri;.!  1 ♦ 1 _ 

6 1 0 C7=D7  09^  A : 

670  IP  \ :9  = 1 THEM  700 
750  IF  1 =0  THEN  MO 
760  IF  R 9 = 1 THEM  790 
780  GOTO  300 

800  PRINT  "PROJECTILE  LENGTH  ! l" METER" 

810  PRINT  "PROJECTILE  D I AMETEF " 5 D ' METER " 

880  PRINT  ‘BASE  AREA"  ! A "METER  SQUARED  " 

340  FRINT  "SPEED  OF  SOUND " ! A l " VELOC I TV  OF  PR0 jElT IlE  ' * V HE  ER 
350  PRINT  "MACH  NUMBER " ! M" ALTITUDE  ' 5 V " METER  " 

860  PRINT  "VISCOSITY" ! M7 "NEW TON-SEC . METER 4 L" 

390  PRINT  "TRANSITION  REYNOLD’S  NUMBER" j.F 1 

8 1 0 PRINT  "BODY  SURFACE  AREA"  METER  M 

820  PRINT  "PROJECTILE  SURFACE  AREA 3 T " METER r _ ' 

9 3 0 R PINT  ' ITEF.AT 1 0 \ \ C 0 N 3 TANT"  5 D 1 

940  PRINT  "TRANSITION  LENGTH’  I *•>  IS"  ? '.'."METER" 

950  PRINT  "SURFACE  AREA  WITH  LAMINAR  FLOW  Is  ' 5 3 ‘ME 'EM- 

8 6Q  PRINT  "SURFACE  AREA  WITH  TURBULENT  FLO U IS"! 33-3 

.70  PRINT  "REYNOLD'S  NUMBER  BASED  ON  T 0 T A L PRO JEC  TILE  LENGTH 

980  PRINT  "LAMINAR  SKIN  FRICTION  COEFFFIC I ENT  IS" !C1 

890  PRINT  "TURBULENT  Skill  FRICTION  COEFFICIENT  i3";Ul 

1000  PRINT  "ITERATION  CRITERIA  IS"! 05 

1010  PRINT  "DYNAMIC  PRESSURE  I S " ! 09 "NEKTON  METER  T L " 

1020  PRINT  "DENSITY  OF  FREESTREAM  IS" ! Z" KILOGRAM-  ME  TER f J" 

1030  PRINT 

1040  PRINT  "SKIN  FRICTION  DRAG  IS"!  D7 "NEWTONS" 

1050  PRINT  "SKIN  FRICTION  DRAG  COEFFICIENT  IS"!C, 

1060  PRINT  ""  

1170  PRINT  IXXXWXHmXXmxmxmmXAAAA.sftArt.v.^AAA...- 

1 1 SO  STOP 
1190  END 


-El  . 


TABLE  IV 


HP  9830  PROGRAM  TO  CALCULATE  BODY  DRAG;  SAMPLE  OUTPUT 


THIS  PROGRAM  •' ALCALA  Es  PRO  JE»!  T ILL  PUD ^ !l 
THIS  PROGRAM  I A STORED  HI  FILE  0. 


F 0 R •£. IJPEF  S 0 N I C V E L U C I TO  I H 3 F. F T y . 


ALTITUDE  IS  0 


METER: 


PROJECTILE  LENGTH  l.^T  ME  TER 
PROJECTILE  DIAMETER  U. 12593  _ _ ME  TER ^ 

B ft *=; E AREA  U.0124G502"  METER  SUUARED 
C peed  Ij F S 01 J H D 3 4 0 1 E I. 0 C I TV  0 F P P 0 JE'.TILE  1 0 L 0 
MACH*"  NUMBER  3 ALTITUDE  0 METER  __ 

I SOOS I TV  1 . S0OOGE-O5  NEW Tui  I--SE1  . ■ ME  i ER  1 - 
T R A N S I T I 0 N R E V N U L D ' S \ < U M BER  2 O 0 0 0 0 0 _ 

BODY  SURFACE  AR  EA  O . E 1 32 1 S 1 3 j ME1 ER  r 2 
P R u J ECTILE  S U P F A C E AREA  0 . S 1 3 2 i U 1 3 9 M E T E P - 

ITERATION  CONSTANT  100  _ MCTCr, 

TRANSITION  LENGTH  > L * » IS  U. 023*1 152b  METcR 

SURFACE  AREA  NITh  LAllINAP  FLOW  lb  0.01  l4‘.'L.-"5  :•  _ L 
SURFACE  AREA  NITH  TURBULENT  FLOW  IS  O.SUlSlb-L* 
REYNOLD’S  NUMBER  BASED  ON  TOTAL  rRuJElTILc.  i-cilGTH 
LAM  I NAR  SK IN  FR I CT I ON  COEFFF I C I ENT  I S _ S . 2c 4 1 2E-04 _ 
TURBULENT  SKIN  FRICTION  COEFFICIENT  Is  1.2490, E U: 
ITERATION  CRITERIA  IS  3. S72B1E-07 
DYNAMIC  PRESSURE  IS  <537245  NEWTON  METER  ^ _ 

DENSITY  OF  FREESTREAM  IS  1.225  K I LUbRhM  METER t :• 


SKIN  FRICTION  DRAG  IS  4S3. 5749443  NEHTUMs 
SKIN  FRICTION  DRAG  COEFFICIENT  IS  O.0BQ37S545 


,:S  .:s . .x.  • 


2.  Tail  fin  Friction  Drag  Calculating  the  jrag  due  to 
friction  on  the  fins  of  Configuration  I was  very  similar  to 
the  methodology  in  calculating  the  body  drag  due  to  friction 
Figure  9 depicts  the  actual  geometry.  Table  J provides  the 
pertinent  dimensions. 


Figure  9.  Tailfin  geometry  for  Configuration  I. 

From  the  data,  the  surface  area  of  each  fin  was  calculated 

A = 2bc  = 0.012097  m2 
s 

The  Reynolds  number  was  calculated  using  the  chord,  c , the 
distance  between  the  leading  and  trailing  edges  measured  in 
the  middle  of  the  tailfin  as  the  characteristic  length 

Re  M -2JL2.  = 4.41  X 10s 

Since  the  Reynolds  number  is  greater  than  2 x 106,  the  fin 
has  both  turbulent  and  laminar  flow.  Calculations  were  made 
to  determine  the  coefficients  of  friction  in  both  the  turbu- 
lent and  laminar  regions.  The  location  of  transition  was 
determined,  and  finally,  the  drag  coefficient  for  each  fin 

was  calculated. 
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using  equations  (12)  through  (21) , the  incompressible 
ana  compressible  shin  friction  coefficients  were  evaluated 
for  the  fins.  The  results  are 


q . = 0.003375 
f i 


and  Cfc  = 0.002056 


in  the  turbulent  region,  and 


Cfi  = 0.000939 


and  Cfc  = 0.000626 


in  the  laminar  region. 

The  location  of  transition  was  found  to  be  at  0.453 

chord  by  using  the  ratio  of  Re*  to  Re. 

The  surface  area  exposed  to  turbulent  and  laminar 
boundary  layer  flow  conditions  was  determined  using  values 

of  chord  and  span. 

A(x*l  * (0.453) (2) (c) (b)  = 0.00548  m2 
ML)  - A (x*)  = (0.547)  (.0121)  = 0.00662  m2 
Drag  was  calculated.  Drag  due  to  the  laminar  flow  region 
was  2.18  N,  and  drag  due  to  turbulent  flow  region  was  8.67  S. 
The  total  drag  per  fin  was  10.35  N.  Finally,  the  drag 
coefficient  was  determined  to  be  0.0014  per  fin  based  on 

fin  area. 

The  above  steps  also  were  programmed  for  the  HR  5830 
computer.  The  program  is  listed  in  Table  V with  a sample 
program  output  in  Table  VI.  By  entering  the  velocity, 
span,  chord  and  altitude,  the  shin  friction  drag  can  be 


calculated. 


TABLE  V 


HP  9830  PROGRAM  TO  CALCULATE  WING  FRICTION  DRAG;  PROGRAM  STEPS 

10  PRINT  "THIS  PROGRAM  CAL.  lULNTE  3 HI  IIC  FRICTION  DPmG. 

HO  PRINT  "THI I PROGPHN  IS  STORED  IN  FILE  2 
30  PRINT  ' " 

40  DEG 
50  D 1=1 00 

6u  DISP  'PRINT  IHF'IJT?YFS=  l NOi  = U"; 

70  INPUT  K 
30  PRINT 

SO  PRINT  FOR  SUPERSONIC  V EL>jl  i TV  INSERT  0.  • 

100  PRINT 

110  DISP  "SEE  PRINTED  STATEMENT  C UNCERnlNG  '</£'■  NT  j ~ ; 

120  INPUT  K9 


1 

30 

D1 

= 

100 

140 

L = 

0 

. 06 

1 

50 

B = 

0 

. OS 

1 

60 

H = 

L 

~B 

1 

70 

Hi 

= 

340 

1 

30 

1.  (l  __ 

if! 

30 

1 

SO 

H = 

i i 

/HI 

00 

H= 

t'' 

620 

2 

10 

i,,i 

1 - 

•*l 

ooo 

22U  PRINT  "■tT+#  + +T*  + T + t**+-»*T  + ft 
230  PRINT  "ALTITUDE  IS " 5 'i'"  MET  PR" 
240  R0= 1.225 
250  Z=RO+ENp'.  -V  H> 

260  M7=l . 3E-05 
270  R 1 =2E+Q6 

2 3 0 =<R1*M?*E F' ■■ : V - H ■ ; ■ « p 0 - V •> 

232  IF  K3=0  THEN  310 

310  S3=2*L+B 

320  IF  L > X THEN  350 

330  3=33 

340  GOTO  420 

350  S=2+6*X 

420  C 1 = 1 . 323  / '-SOROS  l , > 

425  C 1 =C  1 < < 1 +0 . 45*M  1 2 ;•  TO . 25 ) 

4 S W R2  = R 1 * ( L ' X > 

440  PRINT  " 

450  01=C1 

460  Q2=LGT < 0 1 +P2 ) 

470  03=02T2 
4 3 O 0 4=0.053564 ■ 0 3 
4 SO  05=hBS ( 04-0 1 > 

500  IF  Q5<  IE-06  THEN  5F:Fi 
510  IF  05<1 E-05  THEN  570 
520  IF  Q5< IE-04  THEN  550 
530  0 1 =0 1 +C 1 D 1 
540  GOTO  460 
550  01=01+01-"':.  5* D1  > 


TABLE  V (CONTINUED) 


560  GOTO  460 

5 7 tj  Q 1=01+01  40+01.' 

5 ::: O G 0 TO  46 O 

535  Qi  =Q  l ■ •:  < 1 +-0 . 15  IT  L r 0 - 50  1 
530  Q0= PC-  ! V • , . ■ + • 0.5  +£!  F'1  -’V  • H 1 

6 O 0 0 7 - 0. 3 + O'  1 - S +■  0 1 - ' 8 7 - S 1 1 
6 1 0 C 7 - D 7 '■  0!  3 - A ■' 

670  IF  F 1 9=1  THEN  700 
750  IF  K = O THEN  340 
760  IF  K3= 1 THEN  730 
780  GOTO  800 

8O0  PRINT  "WING  CHORD  LENGTH  5 L’TlETER 
810  PRINT  "WING  SPAN" JB" METER" 

880  PRINT  ' BASE  AREA  ' 5 R " ME  TER  3 QUA FED" 

3 4 0 p PINT  " S'  PEED  0 F 3 0 UNOCAL"  V E L 0 C [TV  0 F F PC.  E'  I I L E ? 
850  PRINT  "MACH  NUMBER" ? M "ALTITUDE" 5 V 'METER" 

3 6 O F'R  I N r " V 1 8008  I TV " ^ .17  ' NEWT  ON-SEC . ••  METER  t L ' 

830  PRINT  "TRANSITION  RE',  HOLD’S  NUMBER"?  PI 

310  PRINT  " W IMG  S U R F AC E H R E A " ? S 8 " M E T E P t- L 

980  PRINT  "ITERATION  CONSTANT  '? D 1 

940  PRINT  "TRANSITION  LENGTH" . . + > IS "? 7" METER1 

950  PRINT  "SURFACE  AREA  WITH  LAMINAR  l-  LOW  [ v " < i ' ME  : Ek  * 

960  PRINT  "SURFACE  AREA  WITH  TURBULENT  FLOW  I " * 8 ~ 

970  PRINT  "REYNOLD’S  NUMBER  BASED  ON  TOTAL  WING  CHORD 
980  PRINT  "LAMINAR  SR  IN  FRICTION  COEFFF I C TEN T [S'  ? C l 
9 9 0 PRINT  "TURBULENT  SR  IN  FRICTION  COEFFICIENT  IS"? 01 
1000  PRINT  "ITERATION  CRITERIA  IS"?  0.5 
1010  RR I NT  " DY NAM  I C PRESSURE  I S " ? 09 " NEWT ON.  ME  TER T L " 
ICBC  PRINT  "DENSITY  OF  FREE  STREAM  IS"  ? Z"K  IL0GRAI1- METEh 
1030  PRINT  "" 

1040  PRINT  "SKIN  FRICTION  DRAG  IS" 5 D7" NEWTONS" 

1 050  PR  I NT  " SR  I N FPICT I ON  DRAG  COEFF I C I ENT  I S ' ? 07 
I860  PRINT  "" 

1170  print  " x x x x : •: x x x x x x x x x x x : =: x : < :• : x x x : =: x : i x x : : ; : , x : < : =: x xxxxxx : 

1 1 30  STOP 
1190  END 


MET 


lehc  r 
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TABLE  VI 


HP  9830  PROGRAM  TO  CALCULATE  WING  FRICTION  DRAG 


SAMPLE  OUTPUT 


THIS  PROGRAM  CALCULATES  MING  FRlCTtuN  DRAW. 
THIS  PROGRAM  IS  STORED  IN  FILE 

F f i f SUPERSONIC  VELOCITY  IhS'ERT  O. 


ALTITUDE  IS  3000 


METER 


U IMG  C H 0 R:  D L E N G T H 0 . 0 0 3 5 M E T E R 
WING  SPAN  0.09525  METER 
BASE  AREA  0 . 04333E-0  3 METER  SuUARED 
$ REED  0 F S 0 U N D 3 4 0 , E L 0 C I T < 0 P P R 0 J E C T I L L 
MACH  NUMBER  3 ALTITUDE  3000  METER 
V I S f:  fi  •?.  I T V 1 . S O 0 0 0 E - 0 5 NE  W TON  - S EC.  METER 

TRANSIT  ION  REYNOLD ■ S NUMBER  2000000 
WING  SURFACE  AREA  0.0120907?  METER  f 4 
ITERATION  CONSTANT  100 
TRANSITION  LENGTH?  IT-*  IS  0.042711S44 
SURFACE  AREA  WITH  LAMINAR  FLOW  IS  . L 
SURFACE  AREA  WITH  TURBULENT  FLOW  IS  3 
p E Y N O' L D ' S N U M B E R B A 3 ED  OH  T 0 T A L WINu 
LAMINAR  SKIN  FRICTION  COEFFFICIEHT  !:■ 


1 0 2 0 


METER-  SEi. 


METER 
360  IE-03 
. 3 0 O'  1 4 E — O'  • : 
CHORD  LENGTH 
0 . 204 1 2E-04 


METER T2 


S4i  4.  1 : 


TURBULENT  SKIN  FRICTION  COEFFICIENT  IS  2.1901SE-O3 
ITERATION  CRITERIA  IS  3.57703E-07 
DYNAMIC  PRESSURE  IS  429857.35  NEWTON, METtRtL 
D E N S I T Y 0 F FREE  $ TREA  M I S 0 . 3 4 * S 3 0 3 S 1 K I L U G R h M Mb  i EM  . 


SKIN  FRICTION  DRAG 
SKIN  FRICTION  DRAG 


IS  5.929472104  NEWTONS 
COEFFICIENT  IS  2.2S002E-03 


Tail fin  Wave  Drag 


Wave  drag  was  calculated  for  the  GLM  tailfin.  '~alcu 
lations  were  made  for  double-wedge  airfoils  with  the  aerody- 
namic geometry  as  listed  in  Table  I at  several  angles  of 
attack  and  at  several  Mach  numbers.  Angles  of  attack  between 
0 and  15  degrees  were  analyzed  for  each  airfoil.  Mach  numbers 
considered  were  2.0,  2.5,  and  3.0.  The  procedure  for  calcu- 
lations of  tailfin  drag  was  discussed  in  Section  F of  the 
Aerodynamic  Theory.  Results  are  tabulated  for  each  configura- 
tion in  Section  V. 

4.  Tailfin  Lift  Calculation 

Lift  was  calculated  for  the  GLM  tailfin.  Calculations 
were  made  utilizing  the  double-wedge  airfoil  discussed  in 
Section  G of  the  Aerodynamic  Theory  using  the  aerodynamic 
geometry  listed  in  Table  I for  each  fin.  Methodology  was 
similar  to  the  approach  fcr  calculating  wave  drag  for  each 
fin.  Equation  (29)  is  the  applicable  equation  discussed  in 
the  theory.  Calculations  were  made  for  each  fin  at  angles 
of  attack  between  5 and  15  degrees.  Mach  numbers  considered 
were  2.0,  2.5,  and  3.0.  Results  for  all  three  configurations 
are  included  in  the  lift  coefficients  presented  in  Tables  VIII, 

X,  and  XII  in  Section  V. 

5 . Tailfin  Placement  and  Packaging 

Simultaneously  with  the  calculations  of  fin  drag  and 
fin  lift,  fin  placement  and  packaging  on  the  configuration 
were  analyzed.  In  addition,  thickness  ratio  restrictions, 
fin  deployment  and  storage  were  addressed. 


Since  the  GLM  configurations  are  fired  from  the  5"/54 
MK  45  gun  system,  the  problems  of  packaging  and  storage 
become  apparent.  A location  on  the  body  was  needed  to  allow 
the  fins  to  remain  retracted  during  the  gun  firing  and  then  be 
deployed  upon  exit  from  the  muzzle. 

In  addition,  since  from  the  basic  ground  rules,  the 
propulsion  system  was  to  be  a ramjet,  the  geometry  of  the 
after  body  section  of  the  missile  was  to  be  determined 
primarily  by  the  size  of  the  exhaust  nozzle  required  for 
the  ramjet  to  sustain  Mach  3.0.  Consequently,  fin  placement 
was  bounded  by  the  5"/54  gun  tube  externally  and  the  required 

exhaust  nozzle  size  internally. 

Several  other  parameters  influenced  final  fin  placement 
on  Configuration  I (and  not  coincidentally,  the  subsequent 
configurations) . 

Historically,  the  fin  thickness  ratio,  that  is  the  thick- 
ness to  chord  ratio,  has  been  approximately  5 percent  for  most 
missile  designs.  For  the  tailfins  of  the  GLM  the  wedge  angle 
was  5°.  The  thickness  ratio  was  calculated  to  be  over  8 
percent.  Because  of  the  axial  g-loading  placed  on  the  missile 
by  the  gun  launch,  this  ratio  was  considered  satisfactory  and 
was  not  chanaed.  T this  thesis,  a simple  fin  system  was 
envisioned  for  the  GLM.  The  simplicity  was  required  for  both 
the  packagin  g o£  the  tailfin  system  and  the  mechanism  f 
tailfin  deployment.  Future  study  of  the  tailfin  system  ia 
recommended. 
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Figures  10  and  1.1  show  the  location  chosen  for  tailfin 
placement  on  the  GLM.  Figure  12  depicts  the  envisioned 
design  of  the  hinge  mechanism  used  to  deploy  the  fms  upon 
launch.  The  figures  should  be  viewed  as  sketches  only. 

For  example,  the  inward  tapering  of  the  lines  of  the  GLM 
aft  body  section  in  Figure  10  are  envisioned  only  as  a 
possible  aid  to  keep  the  center  of  pressure  aft  during 

flight. 
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0.09525 


6 . Center  of  Pressure  and  Lif t-to-Draq  Ratio 

Garnell  and  East  [15]  state  that  when  considering 
the  lifting  forces  on  missiles  it  is  convenient  first  of  all 
to  consider  the  combined  normal  forces  due  to  incidence  on 
the  body  and  lifting  surfaces  as  acting  through  a point  on 
the  body  called  the  center  of  pressure.  The  center  of  pres- 
sure location  when  compared  with  the  center  of  gravity  location 
can  be  a good  indicator  of  the  stability  of  a flight  vehicle. 
Calculation  of  the  center  of  pressure  c.p.  location  was 

conducted  for  the  5"/54  GLM. 

In  the  calculation  for  the  location  of  the  center  of 
pressure,  the  lift  coefficient  values  for  the  nose  (including 
the  contribution  by  the  necking) , the  body,  and  the  tailfins 
were  positioned  on  the  configuration  axis.  Figure  13  illus- 
trates the  technique. 


K 


C 


Lt 


GIM  Axis 


Figure  13.  Location  of  lift  coefficients  on  GLM  axis. 
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The  c.p.  location  was  found  by  dividing  the  summation  of 
the  moment  arms  generated  by  individual  components  of  the 
Uft  coefficient  by  the  station  of  the  individual  coeffr- 

cients  of  lift. 


cp 


*-Crn  + Vl,b  * £tCLt_ 

< CLn  + CLb  + CLt  1 


(35) 


, • /tcn  f — 0.6335  in 
After  inserting  values  into  equation  (35) , tcp 

for  Configuration  I. 

Lif t-to-drag  ratio  is  simply,  L/D,  the  lift  divided  by 

the  drag;  the  ratio  also  can  be  determined  by  CL  / C„. 

T an  ancle  of  attack  equal  to  5 , the 
For  Configuration  I at  an  ang 

value  of  L/D  is  5.14  at  Mach  3.0. 

Figure  14  in  Section  V,  shows  Configuration  I;  which  rs 

the  result  of  the  analysis  conducted  in  this  section.  In  a 

similar  fashion,  the  designs  of  Configurations  XI  and  III 

were  calculated.  The  results  also  are  shown. 
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V . AERODYNAMIC  RESULTS 


A.  CONFIGURATION  I 

The  initial  aerodynamic  design  analyzed  in  this  thesis 
was  Configuration  I.  The  major  emphasis  with.  Configuration  I 
was  that  the  inlet  for  the  ramjet  engine  should  utilize  as  much  of 
the  incoming  volume  of  air  as  possible.  Figure  14  depicts 
the  aerodynamic  configuration  while  Table  711  provides  the 
geometric  characteristics  and  Table  VIII  provides  the  calcu- 
lated aerodynamic  data.  Estimated  lift  and  drag  characteris- 
tics are  functions  of  angle  of  attack,  and  the  lift  coefficient 
versus  the  drag  coefficient  is  provided  in  Figures  15  through 
17.  In  order  to  meet  the  loading  system  constraint  specified 
in  Section  II,  Configuration  I was  reduced  in  available  length 
from  1.55  meters  to  1.29  meters.  Because  of  this  reduction 
in  length,  significant  missile  volume  was  lost. 
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TABLE  VII 


GEOMETRIC  CHARACTERISTICS  FOR 


DURATION  I 


BODY : 

Length,  overall,  m. 

Length,  forebody,  m. 
Diameter,  overall,  m. 

Volume,  m3 

2 

Wetted  area,  m 
2 

Base  area,  m 

Cone  angle  of  forebody , deg 


1.3716 

0 

0.1270 

0.0164 

0.5168 

0.0127 

0 


TAILFIN : 

. 2 

0.01210 

Area  (2  panels) , m 

0.0635 

Chord,  m. 

0.09525 

Span,  m. 

Thickness  ratio 

0.0875 

Aspect  ratio 

1.499 

TABLE  VIII 

aerodynamic  data  for  glm  CONFIGURATION  I 


Mach  No. 

2.0 

2.5 

3.0 

(a)  deg 

0.146 

0.122 

0. 105 

0 

D 

n 

0.164 

0.135 

0.115 

5 

D 

O 

0.217 

0.175 

0.143 

10 

D 

C 

0.  305 

0.242 

0.202 

15 

D 

LIFT  COEFFICIENTS: 


Mach  NO. 


2.0 


2.5 


3.0 


-"W 


TABLE  VIII  (CONTINUED) 


AERODYNAMIC  DATA  FOR  GLM  CONFIGURATION  I 


LIFT  COEFFICIENTS: 


Mach  No. 

a = 15 

to 

• 

o 

2.5 

3.0 

% 

0.524 

0.524 

0.524 

% 

1.592 

1.592 

1.592 

% 

0.424 

0.301 

0.280 

=T 

2.540 

2.416 

2.395 

Lif t-to-drag  ratio  at  a = 5 degrees 5.14 

Center  of  pressure,  m. , as  measured  from  nose  tip 

at  Mach  2.0  0.6  34 


B.  CONFIGURATION  II 

Configuration  II  was  the  second  aerodynamic  design 
analyzed  in  this  thesis.  A result  of  the  5 /54  Mk  45  gun 
loading  system  constraint,  the  configuration  was  a compromise 
attempt  to  maintain  the  lowest  possible  nose  drag  on  the 
flight  vehicle  while  maximizing  the  volume  of  the  GLM. 

Figure  18  depicts  the  configuration. 

The  "bottleneck"  forebody  was  designed  to  meet  the 
specific  constraint  of  the  loading  system  to  preclude  any 
system  modification.  The  geometric  characteristics  a.  e 
presented  in  Table  IX,  and  the  aerodynamic  data  are  provided 
in  Table  X.  The  gain  in  overall  volume  from  Conf iguration 
I to  Configuration  II  is  slightly  greater  than  5 percent. 
Lift-to-drag  ratio  was  decreased  from  Configuration  I. 

Lstimated  lift  and  drag  coefficients  as  functions  of 
angle  of  attack  are  presented  in  Figures  19  and  20.  The 
lift  coefficient  versus  the  drag  coefficient  is  depicted  in 
Figure  21. 

During  analysis  of  Configuration  II,  a decision  was 
made  to  modify  the  nose  inlet  configuration.  This  was 
based  on  the  ramjet  propulsion  requirements  in  the  companion 

thesis  of  Brown  [4]. 


TABLE  IX 


GEOMETRIC  CHARACTERISTICS  FOR  GLM  CONFIGURATION 


BODY: 

T^ngth,  overall,  m. 

1.5494 

Length,  forebody,  m. 

0.2889 

Diameter,  overall,  m. 

0.1259 

3 

Volume,  m 

0.0173 

, 2 

Wetted  area,  m 

0.6132 

2 

Base  area,  m 

0.0125 

Cone  angle  of  forebody , deg 

9.5 

TAILFIN : 

2 

Area  (2  panels) , m 

0.0121 

Chord,  m. 

0.0635 

Span,  m. 

0.0952 

Thickness  ratio 

0.0875 

Aspect  ratio 

1.499 

TABLE  X 

AERODYNAMIC  DATA  FOR  GLM  CONFIGURATION  II 


Mach  No. 

2.0 

2.5 

3.0 

(a)  deg 

n 

o 

o 

0.411 

0.330 

0.277 

0 

CD 

0.428 

0 . 343 

0.288 

5 

CD 

0.481 

0.383 

0.320 

10 

cu 

0.569 

0.449 

0.374 

15 

LIFT  COEFFICIENTS: 

Mach  No. 

2.0 

2.5 

3.0 

a - 5 

0.175 

0.175 

0.175 

% 

0.523 

0.473 

0.472 

CL 

lt 

0.164 

0. 150 

0.123 

cT 

0.862 

0.798 

0.770 

0.349 

0.349 

0.349 

0.824 

0.824 

0.824 

0.279 

0.210 

0.183 

0.45  2 

0.383 

1.356 

TABLE  X (CONTINUED) 

AERODYNAMIC  DATA  FOR  GLM  CONFIGURATION  II 


LIFT  COEFFICIENTS: 


Mach  No. 

2.0 

2.5 

3.0 

a = 15 

Cln 

0.524 

0.524 

0.524 

\ 

1.592 

1.592 

1.592 

CL 

0.424 

0.301 

0.280 

lt 

CL 

2.540 

2.417 

2.396 

Lift-to-drag  ratio  at  a 

= 5 degrees 



Center  of  pressure,  m. , 

as  measured 

from  nose 

tip 

at  Mach  2 . 0 


0.839 


C.  CONFIGURATION  III 

Configuration  III  was  the  result  of  the  second  signifi- 
cant compromise  attempted  in  the  thesis.  As  mentioned 
previously,  the  nose  inlet  geometry  of  Configurations  I and 
II  was  not  satisfactory  for  the  ramjet  engine.  Too  much  air 
was  entering  the  combustion  chamber . Configuration  III  was 
the  result  of  the  modification  to  the  nose  inlet  to  reduce 
air  intake.  Figure  22  depicts  the  aerodynamic  configuration 
in  its  entirety.  While  ensuring  that  the  propulsion  require- 
ment was  met  satisfactorily,  care  was  taken  to  ensure  pre- 
viously mentioned  external  geometric  constraint?  still  were 
being  met.  Figure  22  depicts  the  Configuration  III  forebody 
and  the  gun  loading  system  constraint.  Also  depicted  in 
Figure  23  is  the  estimated  volume  for  the  nose  cap  necessary 
to  protect  the  nose  inlet  for  the  missile.  The  resultant 
loss  of  volume  in  Configuration  III  through  modification  of 
the  nose  geometry  of  Configuration  II  was  considered  insig- 
nificant . 

The  geometric  characteristics  for  Configuration  III  are 
presented  in  Table  XI.  Table  XII  provides  the  pertinent 
aerodynamic  data. 

Estimated  lift  and  drag  coefficients  as  functions  of 
angle  of  attack  are  presented  in  Figures  24  and  25.  The 
lift  coefficient  versus  the  drag  coefficient  for  Configura- 
tion III  is  depicted  in  Figure  26.  Lift-to-drag  ratio  was 
decreased  from  Configuration  II  * 
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location  defined  in  Section  IV 


The  center  of  pressure 
was  calculated  for  all  three  configurations.  The  Configura- 
tion III  c.p.  location  was  compared  with  the  center  of 
gravity  location  presented  in  the  companion  thesis  by 
Frazier  [3].  Center  of  pressure  location,  as  measured  from 
the  tip  of  the  GLM  nosev  was  0.0464  meters  aft  of  the  center 
of  gravity  location  measured  from  the  same  point.  Both 
calculations  assumed  GLM  launch  at  Mach  2.0. 
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Figure  22.  Configuration  III. 


TABLE  XI 


GEOMETRIC  CHARACTERISTICS  FOR  GLM  CONFIGURATION  III 


BODY : 


Length,  overall,  m. 

1.5494 

Length,  forebody,  m. 

0.2889 

Diameter,  overall,  m. 

0.1259 

Volume,  m^ 

0.01602 

2 

Wetted  area,  m 

0.6132 

2 

Base  area,  m 

0.0125 

Cone  angle  of  forebody,  deg 

9.5 

TAILFIN : 

2 

Area  (2  panels) , m 

0.0121 

Chord,  m. 

0.0635 

Span,  m. 

0.0952 

Thickness  ratio 

0.0875 

Aspect  ratio 

0.499 

TABLE  XII 


AERODYNAMIC  DATA  FOR  GLM  CONFIGURATION  III 


Mach  No. 

ro 

• [ 

o 

2.  5 

3.0 

(a)  deg 

C°o 

0.529 

0.422 

0.349 

0 

CD 

0.547 

0. 436 

0.  361 

5 

CD 

0.599 

0.477 

0.393 

10 

c 

0.687 

0.542 

0.447 

15 

LIFT  COEFFICIENTS: 


Mach.  No 


2.0 


2.5 


3.0 


TABLE  XII 


AERODYNAMIC  DATA  FOR  GLM  CONFIGURATION  III 


LIFT  COEFFICIENTS: 


Mach  No. 

2.0 

2.5 

3.0 

a = 15 

c 

0.524 

0.524 

0.524 

a 

1.592 

1.592 

1.592 

iJB 

r 

0.424 

0.301 

0.280 

LT 

CL 

2.540 

2.417 

2.359 

Lif t-to-drag  ratio  at  a = 5 degrees 

Center  of  pressure,  m. , as  measured  from  nose  tip 


at  Mach  2 . 0 


0.025 


VX.  reaction  jet  control  theory 


A.  BACKGROUND 

Reaction  jet  control  technology  received  significant 
impetus  in  1964  when,  in  a series  of  wind  tunnel  experiments, 
Zukoski  and  Spaid  [16]  conducted  a study  of  the  flow  field 
around  an  injection  port  for  secondary  injection  of  a gas 
normal  to  the  supersonic  stream.  From  their  data,  they  were 
able  to  calculate  a scale  parameter  that  correlated  the  data 
quite  well.  In  their  analytical  model,  Zukoski  and  Spaid 
likened  the  resultant  shock  pattern  from  the  secondary  infec- 
tion to  that  produced  by  blunt  axisymmetric  bodies.  For  the 
experiment,  the  underexpanded  jet  flow  was  sonic  at  the 
injector,  with  additional  expansion  through  a strong  Prandtl- 
Meyer  expansion  fan.  Wall  boundaries  restricted  the  expansion 
of  the  flow  out  of  the  circular  orifice.  In  1968,  the  data 
bank  was  expanded  by  Zukoski  and  Spaid  [17]  when  they  investi- 
gated the  interaction  of  the  jet  and  freestream  flows  from 
transverse  slots. 

Within  the  next  few  years,  others  [18,  19,  20]  followed 
in  the  investigations  altering  the  parameters  believed  to 

have  an  effect  on  the  freestream  flow. 

In  early  1970,  Wu  [21]  presented  a two-dimensional 

analytical  model  of  the  resulting  force  normal  to  the 
freestream  flow.  The  jet  was  considered  to  be  equivalent  to 
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a solid  body  interrupting  the  freestream  flow. 

Later  that  year,  Nunn  [22]  presented  a method  for  esti- 
mating the  extent  and  shape  of  the  three-dimensional  flow 
field  for  a sonic  jet  exhausting  normal  to  the  supersonic 
freestream.  The  study  was  the  first  to  address  the  problem 
of  jet  interaction  on  bodies  of  revolution. 

Shortly  thereafter,  Grunnet  [23]  found  that  existing 
data  suggested  that  effective  blockage  of  the  equivalent 
flow  body  (two-dimensional)  was  proportional  to  the  inject- 
ant  total  pressure  times  the  area  of  injection  port  flow 
regardless  of  the  injectant  or  freestream  gas  properti.es. 

In  addition,  regardless  of  the  given  freestream  conditions 
and  tiie  injector  port  shape,  the  relative  interaction  force 
was  higher  when  the  gas  was  injected  at  an  angle  upstream  in 
the  freestream  flow.  The  total  side  force  reached  a maximum 
with  the  jet  angled  about  30  degrees  upstream. 

The  interest  in  jet  reaction  control  was  increased  since 
the  early  experiments  seemed  to  confirm  the  belief  that 
under  certain  conditions  the  side  force  produced  by  the 
control  jet  was  greater  than  the  thrust  of  the  jet  alone 
by  some  "amplification"  factor.  Accurate  determination  c£ 
the  amplification  of  the  jet  control  force  resulting  from 
aerodynamic  interactions  between  the  boundary  layer  and  the 
flow  field  over  the  vehicle  associated  with  the  control  jet 
is  the  problem.  Inaccuracy  is  an*  to  the  presence  of  non- 
linear phenomena  associated  with  the  techanics  of  the  inter- 
action; notably  boundary  layer  separation,  viscous  mixing, 
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the  shock-wave  Ooundary  layer  interaction.,  and  a "Mach  disk" 
associated  with  undero/,  jets.  Thus,  only  simplified 

models  for  predicting  j^^'oximatjions  of  the  control  forces 
have  been  developed. 

Define  a Cartesian  coordinate  system  with  coordinates 
x,y,z.  The  wall  lies  in  the  x-z  plane.  The  freestream  flow 
is  in  the  x-direction,  and  the  jet  flow  has  velocity  mainly 
in  the  y-direction.  The  nozzle  is  a slot  which  extends  in 
the  z-direction.  In  a study  conducted  for  the  Navy  by 
Levine  [24],  the  two-dimensional  jet  interaction  flows  using 
the  coordinate  system  discussed  above  were  evaluated.  Scaling 
parameters;  were  developed.  In  contrast  to  three-dimensional 
jet  interaction  such  as  would  occur  with  a circular  nozzle, 
the  high  pressure  zones  in  the  vicinity  of  the  slot  nozzle 
are  constrained  laterally.  The  pressure  distribution  on  the 
wall  is  not  a function  of  z.  In  the  case  of  the  GLM,  the 
z-direction  corresponds  to  a circumferential  dimension  of  the 
body.  The  reaction  jet  control  system  envisioned  for  the  5 /54 
GLM  was  based  on  two-dimensional  jet  reaction  flow  geometry. 

B.  DESIGN  THEORY 

For  the  preliminary  design  of  the  reaction  jet  control 
system  of  the  5"/54  GLM,  the  theory  described  in  the  Levine 
study  [24]  was  used.  The  following  assumptions  and  ground 
rules  applicable  to  the  reaction  jet  were  considered  to  be 
integral  to  that  study. 
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1.  In  order  to  simulate  two-dimensional  flow  mentioned 
previously  and  thus  achieve  the  greatest  amplification  of 
the  thrust  into  a resultant  side  force,  the  injection  slots 
for  the  reaction  jets  should  be  interdigitated  between  the 
GLM  fins 

2.  The  jet  is  underexpanded  and  slightly  supersonic. 

3.  The  reaction  jet  control  slots  are  angled  30  degrees 
upstream  in  the  freestream  flow. 

4.  The  "amplification"  factor  K was  defined: 


where  F is  the  total  measured  normal  force  due  to  the  jet 
thrust  plus  the  aerodynamic  interaction,  and  Fj  is  the 
calculated  jet  thrust. 

5.  The  characteristics  of  the  two-dimensional  jet  inter- 
action flow  field  are  depicted  in  Figure  27  with  an  indica- 
tion of  the  associated  induced  pressure  field. 
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VII. 


CONT ROL  DESIGN  CALCULATIONS 


A.  INTRODUCTION 

The  design  calculations  for  the  hypothesized  reaction 
jet  control  system  in  the  5"/54  GLM  were  divided  into  several 
major  steps.  Utilizing  the  final  iteration  of  the  aerodynamic 
configurations,  -he  coefficient  of  lift  required  for  the 
vehicle  to  execute  a 30 -g  maneuver  was  determined.  From  the 
aerodynamic  results  for  Configuration  III,  an  angle  of  attack 
corresponding  to  the  required  lift  coefficient  was  calculated. 
The  next  step  involved  calculating  the  moments  resulting  from 
each  of  the  individual  lift  coefficients  corresponding  to  the 
nose,  body  and  tail  of  the  GLM.  The  center  of  gravity  loca- 
tion needed  for  the  moment  calculations  was  assumed.  Accurate 
location  of  the  center  of  gravity  is  left  to  future  studies 
of  the  GLM,  From  the  moment  calculations,  the  required  reac- 
tion jet  control  force  was  determined,  and  a graph  depicting 
the  reaction  jet  force  as  a function  of  center  of  gravity 
location  was  constructed.  Finally,  using  the  two-dimensional 
theory  discussed  in  Section  VI,  the  required  F.  (jet  thrust) 
was  calculated.  The  calculated  jet  thrust  is  the  value 
corresponding  to  the  mass  flow  bled  off  the  ramjet  engine  for 
GLM  control.  The  reaction  jet  control  design  is  sketched  in 

Figure  31. 
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B.  REQUIRED  LIFT  COEFFICIENT 

One  of  the  design  specifications  for  the  GLM  was  that 
the  missile  be  capable  of  a 30 -g  maneuver  to  intercept  a 
target.  The  first  step  in  the  calculation  of  the  required 
reaction  jet  force  needed  to  maneuver  the  GLM  is  to  determine 
the  required  lift  (or  normal)  force  coefficient  for  the 
missile.  As  shown  in  the  Appendix,  lift  is  nearly  equal  to 
the  normal  force  N when  the  angle  of  attack  is  small  and 
the  lift-to-drag  ratio  is  large.  Since  the  conditions  of 
small  a and  large  L/D  apply  to  the  GLM,  L and  N are  used 

interchangeably  in  this  thesis. 

From  Figure  28  the  lift,  L,  required  by  the  GLM  in  a 
constant  acceleration  turn  is 

L = nmg 

The  maneuver  load  factor  is  n,  and  m is  the  missile  mass,  kg. 


L 


Figure  28.  GLM  in  a constant  acceleration  turn 


Recall  that 


L = CL  qSr  (35) 

S is  the  reference  area  for  the  GLM.  The  required  lift 
r 

coefficient  is 

c „ nrn9  (36) 

Lr  q Sr 

Equation  (36)  gives  the  lift  coefficient  needed  to  develop 
the  required  lift  for  the  GLM  to  maneuver.  Since  the  GLM  is 
considered  to  be  powered  to  intercept,  complete  burn  of  the 
rocket  fuel  and  0.85  burn  of  the  ramjet  fuel  were  used  in  the 
calculation . For  the  calculations,  values  for  the  mass  of 
rocket  fuel  and  the  mass  of  the  ramjet  fuel  were  obtained 
from.  Brown's  [4]  investigation  of  ramjet  propulsion.  The 
weight  of  the  rocket  exhaust  nozzle  was  used  in  the  calcula- 
tion since  upon  complete  burnout  of  the  rocket  fuel  the 

nozzle  is  expended. 

Mass  (Ramjet  fuel)  = 3.326  kg 
Mass  (Rocket  fuel)  - 5.2512  kg 

Mass  (Nozzle)  = 0.5  kg 

Firing  mass  of  the  GLM  is  47.14  kg.  By  subtracting  the 

values  calculated,  the  maneuver  mass  is  38.56  kg. 

Design  values  for  maneuver  at  sea  level  were  inserted 
into  equation  (36).  calculated  to  be  1.427. 

From  the  lift  curve  shown  in  Section  V for  Configuration 
XU,  Figure  19,  an  angle  of  attack  of  11  degrees  is  required. 


c.  CALCULATION  OF  MOMENTS 

The  second  step  in  the  calculation  of  the  required 
reaction  jet  force  involved  the  calculation  of  the  moments 
coincident  to  each  of  the  individual  lift  coefficients 
corresponding  to  the  nose,  body  and  tail  of  the  GLM. 

Figure  29  shows  the  location  of  the  individual  component 
lift  coefficients  for  Configuration  III.  Distances  m the 
calculations  were  measured  from  the  nose  tip  of  the  oLM. 


Figure  29 


Lift  (normal  force)  coefficient  locations  for 
Configuration  III. 


Because  the  reaction  jet  force  was  to  be  interdigitc'ced 
between  the  tailfins,  the  distance  of  this  fore*  fvom  the 
nose  tip  was  equated  to  that  of  the  component  f'.'fCC*  due 
to  the  tailfins.  The  component  of  lUt  due  to  the  co.nf rea- 
ction "bottleneck"  was  incorporated  i*  to  the  cr*  —nr.it  of 
force  due  to  the  nose  itself,  as  discuss*  clie  ''htc 


The  overall  length  of  the  GLM,  L,  and  measured  distances 

as  depicted  in  Figure  29  from  the  nose  are: 

L = 1.5494  m. 

= 0.1016  m. 

I = 0.7049  m. 

B 

? = 1.5176  m. 

T 

The  center  of  gravity  location  xcg  was  assumed  to  be 

one  half  the  length  of  the  GLM  for  the  initial  calculation. 

For  more  accurate  determinations,  a plot  of  the  force  due 

to  the  reaction  jet  required  to  trim  the  GLM  for  a 30-g 

maneuver  at  Mach  3.0  and  at  sea  level  was  made. 

x = 0.7747  m. 
eg 

From  the  Configuration  III  aerodynamic  data,  the  compo- 
nent force  coefficients  for  lift  were  determined  to  be 

C = 0.384 

h 

c_  = 0.978 

LB 

CT  = 0.299 

LT 

for  the  required  angle  of  attack  a for  the  30-g  maneuver. 

For  each  component  force  coefficient  of  lift,  the 

corresponding  moment  coefficient  CM  was 

C . = CT  • ( xra  - l,  ) (37 

cMj  10  eg  3 

where  the  index  j indicates  nose,  body,  tail  or  reaction 
jet.  inserting  the  values  for  each  component  of  lift  into 

equation  {3n)  : 
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0.2585 


C,  = 0.0683 


"Mr 


= -0.2221 


For  the  GLM  to  be  trimmed  at  the  angle  of  attack  required 
for  the  30 -g  maneuver 


+ % + 


+ C “ ° 

T CMrj 


(38) 


where  C 


'M 


is  the  moment  coefficient  required  of  the  reaction 


ro 


jet  side  force. 

By  retracing  the  calculations  for  the  given  location  of 
the  reaction  jet  control  orifices,  the  side  force  can  be 
found.  Using  equation  (37) 


and 


cr  = CN  ” °*141° 

Lrj  Nrj 

Force  . = C qSr  = 1122.6  N 
rJ  r j 


For  a given  center  of  gravity  and  reaction  jet  moment 
coefficient  required,  the  reaction  jet  force  has  been  calcu- 
lated.  The  required  forces  due  to  the  reaction  jet  have 

been  plotted  in  Figure  30. 
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D.  CALCULATION  OF  REACTION  JET  THRUST 

From  the  theory,  recall  that  the  side  force  normal  to 
the  bcdy  was  not  equal  to  the  actual  reaction  jet  force  but 
in  fact  was  greater  by  an  amplification  factor  K [21,  24]. 
Based  on  the  tabulated  results  from  the  Levine  study  [24], 
the  amplification  factor  for  the  GLM  design  was  conservatively 
assumed  to  be  two.  Nunn  [22]  states,  however,  that  three  as 
an  amplification  is  not  unreasonable  for  the  two-dimensional 
approximation.  The  reader  is  referred  to  Refs.  [21-24]  for 
more  information.  Thus,  for  K = 2 the  reaction  jet  thrust 
needed  for  the  30-g  maneuver  is 

F.  = F ./K  = 561.3  Newtons 

3 ry 

The  jet  thrust  is  that  force  which  is  result  of  the  mass  bled 
from  the  ramjet  engine.  Figure  31  depicts  a preliminary 
design  sketch  of  the  reaction  jet  control  orifices  interdigi- 
tated  between  the  tailfins.  Orifices  would  be  located  45 
degrees  between  fins. 


Reaction  jet  control  orifices  located  between 
tailfins  of  Configuration  III. 
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Figure  31. 


VIII.  CONCLUSIONS 

A,  AERODYNAMIC  CONCLUSIONS 

The  following  aerodynamic  conclusions  were  reached  as 

a result  of  this  study: 

1.  The  drag  coefficient  increased  from  Configuration  I 
through  Configuration  II  to  Configuration  III.  Lift-to- 
drag  ratio  decreased  from  Configuration  I through  II  to  III. 

2.  Gun-launched  missile  volume  increased  from  Conf iguratioi 
I to  II  but  decreased  slightly  in  Configuration  III . 

3.  The  calculated  center  of  pressure  was  found  to  be  aft  of 
the  estimated  center  of  gravity  for  a Mach  2.0  launch  from 
the  gun.  The  margin  of  stability  was  3 percent  of  the  length 

of  the  GLM. 

B.  REACTION  JET  CONTROL  CONCLUSIONS 

The  following  reaction  jet  control  conclusions  were  reached: 

1.  The  total  force  of  reaction  control  jet  thrust  required 
to  trim  the  GLM  to  11  degrees  angle  of  attack  for  a 30-g 
maneuver  at  sea  level  Mach  3.0  is  561.3  Newtons.  The  amplifi- 
cation was  assumed  to  be  2 . 

2.  The  static  pressure  in  the  chamber  for  Mach  3.0  as 
determined  by  Brown  [4]  is  2,137,873  N/m  . 

3.  The  percentage  bleed  of  the  ramjet  airflow  required  for 
the  reaction  jet  control  is  19  percent.  Ramj  t design  did 
not  include  this  bleed. 


96 


APPENDIX 


[ 


RELATION  BETWEEN  AERODYNAMIC  FORCES 

The  purpose  of  this  appendix  is  to  clarify  the  relation 
between  aerodynamic  forces  on  a flight  vehicle.  Figure  32 
helps  depict  the  situation. 


t4 


Figure  32.  Relation  between  aerodynamic  forces. 


In  Figure  32,  the  following  symbols  are  defined: 

V - freestream  velocity  vector 
00 

L - aerodynamic  lift  which  is  normal 
to 

D - aerodynamic  drag  which  is  parallel 

to  V ' 

00 

N - aerodynamic  normal  force  which  is 
normal  to  the  missile  axis 


A - aerodynamic  axial  force  which  is 
parallel  to  the  missile  axis 

T - total  aerodynamic  force 

cl 
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From  the  geometry  of  Figure  32, 


T = L2  + D2  = L A + (D/L)  2 
a 


(A-l) 


and  likewise 

8 = tan"1  D/L  (A_2) 

If,  from  the  geometry,  the  aerodynamic  normal  force,  N,  which 
is  normal  to  the  missile  axis  is  equal  to  the  total  aero- 
dynamic  force  Ta  multiplied  by  the  cosine  of  6 sinus  a then, 
by  substituting  equation  (A-l)  and  by  the  use  of  a trigono- 
metric identity 

N = L /l  + (D/L)  2 (cosBcosa  + sinBsina)  (A-3) 

But  from  the  geometry  of  Figure  32,  equations  !A-4)  and  (A-5) 


can  be  seen. 


COS0  = 


( A-4 ) 


XT 


(D/L) 


sin8  = 


D/L 


(A-5) 


XT 


(D/L) 


By  substitution  of  equations  (A-4)  and  (A-5)  into  equation 
(A-3) , equation  (A-6)  can  be  found. 


N =■  L 


XT 


(D/L) 


cosa 


(D/L) sina 


XT  (D/L)2  A + (D/L)2 


and 


N = L (cosa  + D/Lsina) 
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Hence  for  small  a 
axial  force  A is 


and  large  L/D,  N - L. 


From  Figure  32  the 


Combining  equations 


For  situations  of  a 


A = T sin  (S-ct) 
a 

(A-l).  (A- 4) , (A-5)  and  (A-8)  yields 

A = l (D/L  cosa  - sina)  (A“ 

small  angle  of  attack  a and  a large 


L/D,  A » D. 
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